Zfe(10650) and Zf,(10610) states in a chiral quark model 
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We perform a systematic study of BB* , B* B* , DD* and D* D* systems by using effective in- 
teraction in our chiral quark model. Our results show that the interactions of BB* , B* B* , DD* 
and D*D* states are attractive, which consequently result in BB* , B* B* , DD* and D*D* bound 
states. The recent observed exotic-like hadrons of Zb(10610) and Zi,(10650) are, therefore in our 
approach, interpreted as loosely bound states of BB* and B* B* , while X(3872) and y(3940) are 
DD* and D* D* molecule states, respectively. 

PACS numbers: 13.75.Lb, 12.39.-x, 14.40.Rt 



I. INTRODUCTION 



Recently, the Belle Collaboration pronounced two 
charged bottomonium-like resonances Zb(10610) and 
Zb(10650) [l]. These two resonances were observed in the 
five decay channels of T(5S') — > 7r='=T(nS') [n = 1,2,3) 
and T(5S') — > Tr^ht{mP) (m = 1,2). The recommended 
masses and widths of the two resonances are M(10610) — 
10608.4 ± 2.0 MeV/c2, Fi = 15.6 ± 2.5 MeV/c^ and 
M(10650) = 10653.2 ± 1.5 MeV/c^ Ta = 14.4 ± 3.2 
MeV/c^, respectively. For both states, experimental re- 
sults fixed their quantum numbers as /'^( J^)=l^(l^). 
The two resonances can't be simply interpreted as the 
bb heavy quarkoniums because they are charged. It is 
suggested that they might be BB* and B*B* molecules, 
respectively, since their masses are near to the thresholds 
of BB* and B*B* , and they have narrow widths. 

There are many different explanations about the struc- 
tures of Z&(10610)and ^,,(10650) in the literature [3- 
[Ts! ]. Ali et al. 0, [1] and Guo et al. 0] supported 
that Zb(10610) and Zb(10650) could be explained as 
tetraquark states. Cui et al. @ showed both tetraquark 
and molecule interpretations were compatible with the 
experimental measurements in the QCD sum rule cal- 
culation. However, Matheus @ and Chen Q argued 
that the masses of tetraquark states were obviously below 
the BB* and B*B* thresholds. Zhang et al.[8.] showed 
their quite satisfactory result that Zf,(10610) might be re- 
garded as a BB* molecule which agreed with the Belle's 
molecule suggestion Moreover, Nieves and Valder- 
rama argued that there might be three isoscalar BB* 
bound states by considering one pion exchange interac- 
tion, but all of the three didn't fit the quantum numbers 
of Zb(10610) [9^. Ding et al.[l3IIl| suggested BB* and 
B*B* could form molecule bound states, but these states 
couldn't fit to the 1+(1+) quantum numbers of the two 
Zf) states. By including only screened color-Coulomb and 
screened linear confining potentials, Wong [T^] also got 
some BB* and B*B* molecules, but his Zf,(10610) arid 
Zf,(10650) should be the first radial excitation of BB* 
and B*B* molecules, respectively. In addition, by con- 
sidering both one boson exchange interaction and the 
one pion exchange interaction and considering the S ~ D 



mixing between B/B* and B* , Sun et al. [ij] gave their 
result of Zb(10610) and Zb(10650) in the molecule sce- 
nario. It should be mentioned that the calculations of 
Dong et al. [l3| and Cleven et al. [15] supported the 
molecule scenario, however, the latter couldn't exclude 
the tetraquark component. 

We see that by using different approaches one may 
have different interpretations for the two new resonances 
of Zf,(10650) and ^^(10610). In this work, we will study 
the two resonances based on our chiral quark model. It 
should be stressed that the chiral quark model is success- 
ful to well describe the binding energies of the baryon 
ground states, the binding energy of the deuteron, the 
NN scattering phase shifts and the NY cross sections 
simultaneously [16 25]. It is expected that the chiral 
quark model provides another way to study the possible 
BB* and B*B* bound states problem on quark level. 

In addition, it is known that the study of the hadronic 
states of X(3872) observed in Ref. [lill and Y(3940) 
observed in Ref. [33-^3l have been a hot issue recently. 
This is because they have exotic properties and can't be 
simply regarded as cc charmoniums. Among the vari- 
ous approaches, hadronic molecule of 0^(1^"'') DD* is 
one feasible interpretation of X(3872) (sj-Ilil. _There- 
fore, analogous to the systems of BB* and B*B*, it is 
also of significance to study the DD* and D*D* systems 
in our chiral quark model for X(3872) and Y(3940). 

It should be mentioned that several works based on 
our chiral quark model have been published [l6l - [l8| to 
study the bound states and binding energies of heavy- 
meson molecule systems with Resonating Group Method 
(RGM). The RGM method is a successful one proficiently 
adopted in Ref. [T6l - l25| . In this work, different from the 
RGM method, we apply the approach introduced in Ref. 
poj to derive the analytical from of the total interaction 
potentials between the two S-wave heavy mesons. Here 
the interaction potentials result from the one meson ex- 
change interactions between two light quarks. We antic- 
ipate that this method would give a more accurate de- 
scription of the short-range interaction between the two 
clusters in the chiral quark model, and we thoroughly 
investigate the possible bound states of BB* and B*B* 
systems as well as their charm partners DD* and D*D* 
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by solving the Schrodinger Equation with our analytical 
interaction potentials between the two clusters. 

The paper is organized as follows. In section the 
framework of our chiral quark model is briefly introduced, 
and the analytical forms of the effective interaction po- 
tentials between the two S-wave meson clusters under our 
chiral quark model are given. The bound state solutions 
for the BB* and B*B* systems and their charm part- 
ners of DD* and D*D* are shown in Sec. Hill Finally, 
the summary is given in Sec. IIVI 



II. FORMULATION 

The framework of our effective chiral quark model has 
been discussed extensively in the literature [H, [13, IT9I - 
[2ll |. Here we just briefly present the main part of the 
model. The total Hamiltonian of the approach reads 



(1) 



where Tq is the sum of the kinetic energy operator of 
all the quarks and antiquarks in the two clusters, and 
Tg is the kinetic energy operator of the center-of-mass 
motion. In a molecule scenario or cluster model, because 
the interactions inside each cluster don't dedicate to the 
total interaction between the two clusters of the system, 
only the interactions between the two different clusters 
need to be taken into account. V{rij) in Eq. 1 represents 
the interactions between the ith light quark or antiquark 
in the first cluster and the jth light quark or antiquark 
in another one, and 



For light quark-quark interaction, V"'''{rij), V'^''{rij) 
and V'"(rij) have such forms: 

V-inj) = -C(gc/.,m,„,A)Xi(m,„,A,ry)(A°A°),(6) 



TTl 

^'^"(^y) = C:{gch,mj,^,A)-—-^ — X2(m^^,A,ry) 



(7) 



Gmiirij 



5chv Mn g^hv ^^N 



(X'^X^) , (8) 



where A° is the Gell-Mann matrix in flavor space. In our 
chiral SU(3) quark model, we only consider the scalar me- 
son and pseudoscalar meson exchanges, while in our ex- 
tended chiral SU(3) quark model vector meson exchange 
interactions are also included, m^^ , m^^ and mp^ in Eq. 
5 are the masses of the scalar nonets, the pseudoscalar 
nonets and the vector nonets, respectively. Mn in Eq. 
8 is a mass scale usually taken as the mass of nucleon 
[2H . In the above equations, rrii is the ith light con- 
stituent quark mass in the first cluster, rrij is the jth 
light constituent quark mass in the second cluster, gch 
is the coupling constants for the scalar and pseudoscalar 
nonets, gchv and fchv are the coupling constants for the 
vector coupling and tensor coupling of vector nonets, re- 
spectively. In Eqs. milllHl 



Here V^'~^^ is the one-gluon exchange (OGE) potential 
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3 nig^mg^ 



(3) 



and Vff"'^ is the confinement 



^™"^(r-.,) = -(A?-A^)(a^r,,-Ka^) 



(4) 



V (rij ) in Eq. 2 stands for the total one-meson ex- 
change interaction potential. Generally, 

V'ir..,) = ^ V^" in,) + ^ V^' in,) + j2 {n,),{5) 



with y {rij ) , {rij ) and y " (r^ ) being the interac- 
tions induced from scalar meson, pseudoscalar meson, 
and vector meson exchanges, respectively. 



C{gch,m,A) 
Xi{m,A,rij) 
X2{m,A,rij) 
Y{x) 



A2 



2™' 

Y{mr,,)~-Y{Ar,,), 



(9) 
(10) 



Y{mr,,)-[ A ) r(Ar„), (11) 
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For light quark-antiquark interactions, G, 



(12) 



the 



G-parity of the exchanged meson, should be included 
in the interaction potentials V"'^{rij), V^'^{rij) and 
V'>^{r.,). 

In the BB*, B*B*, DD* and D*D* two-heavy-meson 
systems, only the light quark-antiquark interactions are 
considered. To solve the bound state problem for these 
systems, unlike the Resonating Group Method (RGM) 
employed in our previous works Il6|, 1711 . we use the ap- 
proach explicitly discussed in Ref . [40] . We first derive an 
analytical form of the total effective interaction between 
the two mesons, and then by solving the Schrodinger 
Equation, we calculate the binding energy E, and finally 
draw a conclusion whether the molecule bound states 
could exist in the two-meson systems or not. 
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Now corresponding to scalar, pseudoscalar and vector 
meson exchanges (see Eqs. HI [7] and |S]) , using the method 
described in Ref. [13], we get the analytical effective in- 
teraction potentials between the two S-wave heavy meson 
clusters as 



V^iiO = -G,„C(5ch,m,„,A)Xi,,-(m,„,A,O(A^A|0l3) 

9 

m 



(14) 



+ - 



fchv mq + niq ^ /c/i^, mqmq 
eniqiTiq y ' gchv Mn glt^^ Mil 



1 



xX2qq{mp^,K,i) {aq ■ ffq) 



(15) 



Here, ^ is the relative coordinate between two different 
clusters, namely, the relative coordinate between the two 
centers-of-mass of the two clusters, and 

Xlqq{m,A,0 = YqqimO--YqqiA^), (16) 



with gJij^^/An = 13.67 determined from experiments. In 
our extended chiral SU(3) quark model, vector meson 
exchange interactions are included. The coupling con- 
stant for vector coupling QcHv hi Eq. |S] is taken to be 
2.351 and 1.973, respectively. Correspondingly, the ratio 
of the tensor coupling fchv to the vector coupling gchv is 
taken to be and 2/3, respectively. In our calculation, 
the masses of the mesons are taken from the PDG 41], 
except the a meson, which does not have a well-defined 
value. Here nia- is obtained by fitting the binding energy 
of the deuteron [2lj. It is Too- — 595 MeV in our chiral 
SU(3) quark model, 535 MeV and 547 MeV for fchv /gchv 
respectively taken as and 2/3 in the extended chiral 
SU(3) quark model. The cutoff A is taken as 1100 MeV. 
Since there is no color-interrelated interaction between 
the two color-singlet clusters, we don't have to consider 
the OGE and confinement interactions, and therefore, we 
don't list the OGE coupling constant gij in Eq. [3] and 
the parameters of the confinement potential a^j and a^j* 
in Eq. H 

The remaining parameters to be determined are the 
heavy quark masses TOc and mi,. In our work, we find 
that the final results are not sensitive to the variation of 
the heavy quark masses, and we take TOc = 1430 MeV 
nil and TOb = 4720 MeV [Hi as typical values. 



X2qq{m,A,0 = Yqq{mO-{^] YqqiA^- (17) 



III. NUMERICAL SOLUTIONS 



The modified Yukawa term reads 



2to^ 



Here, 



/3 = 



f^Qq 



(18) 



(19) 



and Lo is the harmonic-oscillator frequency of the light 
quark wave function, q (q) and Q (Q) in Eq. 19 are the 
light quark (antiquark) and heavy quark (antiquark)in 
the two clusters, respectively. 

There are some necessary parameters in the interac- 
tion potentials in our chiral quark model. In this work, 
we adopt the parameters determined in our previous 
works [id, [H Elil. The up/down quark mass Too 
is taken as 313 MeV. The harmonic-oscillator frequency 
w is 2.522/to~^ in the chiral SU(3) quark model and 
3.113/to~^ in the extended chiral SU(3) quark model. 
The coupling constant for the scalar and pseudoscalar 
chiral field gch — 2.621 is fixed by the relation of 



To make sure the BB* , B*B* systems and their charm 
partners of DD* , D* D* have definite quantum numbers 
of isospin / and C-parity C, we construct the following 
flavor wave functions as Refs. [H [ll|39|. For BB* 
system: 

( 1 

(B*+ B" + cB+ B*'^) 

r 

I = 1 : < —{B*- B" + cB' B*^) (20) 

r 

-[{B*^B° - B*+B-) + c{B"B*° - B+B*-)], 
< 2 

/ = : ^[{B*°B^ + B*+B-) + c{B°B*^ + B+B*-)]. (21) 



Here, c=l for C=+ and c=-l for C=-. For B*B* sys- 
tem, the flavor function is 



/ = 1 




-^{B*"B*° - B*+B*-), 
v2 

/ = 0: -^(B*^B*° + B*+B*-). 
V2 



(22) 



(23) 



9ch _ 9 9nn^ 



An 25 47r M 2, ' 



For DD* and D*D* systems, the similar expressions can 
be obtained. 
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BB* 



By using the analytical effective interactions we have 
gotten in section II, we study the BB* system with 
different isospin / and C-parity C. After solving the 
Schrodinger Equation with the programs developed in 
Refs. [iJiill and with the interactions in Eq. I13I14I151 
we list the obtained results in Table HI In the table, I, 
II and III refer to the chiral SU(3) quark model, the 
extended chiral SU(3) quark models with vector nonets 
coupling constant fchv/gchv=0 and 2/3, respectively. If 
bound state does exist, we list, in the table, the binding 
energy E in MeV. The symbol " — " denotes that bound 
state doesn't exist. 

TABLE I: Binding energies (MeV) of the bound states of BB* 



I 


c 


I 


II 


III 





+ 


46.74 


97.75 


78.04 




3.27 


66.53 


50.59 


1 


+ 


0.39 


1.75 


1.49 



From Table HI we see there are three 5-wave BB* 
bound states with quantum numbers /'^(J^'") — 
0+(l++), 0-(l+-) and 1+(1+") respectively. The bind- 
ing energies of 0+(l++), 0~(1"' ) and ) bound 
states are respectively 46.76-97.75 MeV, 3.27-66.53 MeV 
and 0.39-1.75 MeV in our models. Table I shows that 
the attractive interaction of 1=0 BB* system is much 
stronger than that of 1=1 BB* system and the attrac- 
tive interaction in the extended chiral SU(3) quark model 
(models II and III) is stronger than that in the chiral 
SU(3) quark model (model I). From our calculation, we 
notice that in BB* system, a, a' , tt, oj, and p exchange 
interactions play a dominant role and they determine 
whether BB* bound state exists or not. 

For the /=0 C—+ case, in the chiral SU(3) quark 
model, (7, cr' and tt exchanges provide strong attraction, 
so the total interaction is strong enough to form a BB* 
bound state. In the extended chiral SU(3) quark model, 
the contributions of vector meson exchange are also in- 
cluded, and p and w exchanges provide additional attrac- 
tion. Therefore, the BB* system has a larger binding 
energy as shown in Table I where this 0+(l++) bound 
state has a binding energy 46.76-97.75 MeV. 

For the 7=0 C=- case, in the chiral SU(3) quark 
model, a and a' exchanges provide strong attraction but 
TT exchange provides strong repulsion, as a result of the 
competition, the BB* still binds but the binding energy 
becomes much smaller. In the extended chiral SU(3) 
quark model, p and w exchanges provide very strong at- 
traction as well, so the binding energy becomes larger 
than the one in the chiral quark model. This 0~(1"' ) 
bound state has a binding energy 3.27-66.53 MeV. 

For the 1=1 C=+ case, in all the three quark mod- 
els, (7 exchange provides strong attraction but a' and tt 



exchanges provide strong repulsion, thus the total inter- 
action is weakly attractive, however the BB* can't bind 
in our approach. 

For the 1=1 C=+ case, in the chiral SU(3) quark 
model, a exchange provides strong attraction but a' and 
TT exchanges provide strong repulsion, thus the total in- 
teraction is weakly attractive, however the BB* can't 
bind in our approach. In the extended chiral SU(3) quark 
model, the strong attraction provided by tu exchange and 
the strong repulsion provided by p exchange almost can- 
cel each other, so the BB* still can't bind. 

For the 7=1 C=- case, in the chiral SU(3) quark 
model, (T and tt exchanges provide strong attraction but 
cr' exchange provides repulsion, so the total interaction is 
attractive, and the BB* can form a weakly bound state. 
In the extended chiral SU(3) quark model, the contri- 
butions of vector meson exchange are also cancelled, so 
the BB* still can form a 1+(1"' ) weakly bound state 
with weak binding energy 0.39-1.75 MeV. In Fig. [H the 
total interaction potential of the charged 1+(1+) S'-wave 
BB* system is depicted. This charged 1+(1+) S'-wave 
BB* weakly bound state has a mass 10603-10604 MeV, 
and it might be explained as one possible candidate of 
Zfc(10610). 




FIG. 1: The total interaction potential of 7=1 and C=- BB* 
system. The solid, dashed and dotted lines represent the re- 
sult obtained from the chiral SU(3) quark model, the extended 
chiral SU(3) quark model with fchv/Ochv taken as and 2/3, 
respectively. 



B. B*B* 

By solving the Schrodinger Equation with the analyti- 
cal potentials of B* B* system with different isospin 7 and 
spin S, we list our results of the bound state problems of 
the system in Table HIl 

According to Table HTJ we see that four S-wave bound 
states with different quantum numbers could exist in 
B*B* system, namely, the bound states of 1+(1^ ), 
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TABLE II: Binding energies (MeV) of the bound states of 
B* B* system with different I and S quantum numbers. Sym- 
bols have the same meaning as in Table HI 



/ 


s 


I 


II 


III 









51.65 


37.55 





1 


3.30 


66.73 


50.77 




2 


46.91 


97.98 


78.24 







1.51 


3.66 


3.32 


1 


1 


0.41 


1.78 


1.52 




2 









with binding energy 0.41-1.78 MeV, 0^(1+") with bind- 
ing energy 3.30-66.73 MeV, 0+(2++) with binding en- 
ergy 46.91-97.98 MeV and 1^(0++) with binding energy 
1.51-3.66 MeV. From our calculation, we notice that, in 
B*B* system, a, a', tt, uj and p exchanges also play a 
dominant role to determine whether B* B* could form a 
bound state or not. 

For the /=0 5=0 case, in the chiral SU(3) quark model 
(model I), a and a' exchanges provide attraction, how- 
ever, TT exchange provides strong repulsion. As a result 
of this competition, the total interaction is too weak to 
bind B*B* . In the extended chiral SU(3) quark model 
(models II and III), because the contributions of vector 
meson exchange are also included, and p and uj exchanges 
also provide very strong attraction, the B*B* can form 
a bound state. 

For the 7=0 S=l case, in the chiral SU(3) quark model, 
a and a' exchanges provide attraction but tt exchange 
provides repulsion, and the total interaction in the B*B* 
system is still strong enough to form a bound state. In 
the extended chiral SU(3) quark model, the additional p 
and UJ exchanges provide very strong attraction, there- 
fore, the binding energy becomes larger. This 0~(1"' ) 
bound state has a binding energy of 3.30-66.73 MeV. 

For the /=0 S—2 case, in the chiral SU(3) quark model, 
all fj, a' and tt exchanges provide strong attraction, so the 
B*B* forms a bound state. Moreover, in the extended 
chiral SU(3) quark model, p and lu exchanges provide ad- 
ditional attraction, so the B*B* system could be deeply 
bound. This 0+(2++) bound state has a binding energy 
46.91-97.98 MeV. 

For the two cases of 7=1 5=0 and 1, in the chiral 
SU(3) quark model, a and tt exchanges provide attrac- 
tion, and a' exchange provides repulsion. However, the 
total interaction is attractive, and the B*B* can form a 
weakly bound state. In the extended chiral SU(3) quark 
model, like the 7=1 BB* system, the contributions of 
vector meson exchange are also almost cancelled, and the 
B*B* still can form weakly bound states. The 1~(0++) 
and 1"'"(1"' ) bound states have weakly binding energies 
1.51-3.66 MeV and 0.41-1.78 MeV, respectively. In Fig. 
m the total interaction potential of the charged 1"'"(1''") 
5- wave B*B* system is plotted. It should be stressed 
that this 1"'"(1''") B*B* weakly bound state has a mass 
10648-10650 MeV, which is consistent with the mass and 



quantum numbers of ^{,(10650). As a result, the new 
resonance of Zb(10650) might be explained as the weakly 
bound state of 1+(1+) B*B* in our approach. 




-30 - 
-35 - 
-40 - 



0.0 0,5 1.0 1.5 2,0 2,5 

r(fm) 

FIG. 2: The total interaction potential of 7=1 S=l B* B' 
system. Symbols have the same meaning as in Fig. [T] 

Comparing Fig. [2] with Fig. [TJ we see the two interac- 
tion potentials are almost the same. It is not surprising 
because in these two cases, the spin matrix elements in 
Eqs. [14] and [15] are both -1 and the other parts of the 
interaction potentials are the same as each other. The 
B*B* system has a slightly larger binding energy because 
of its larger reduced mass than the 737?* system which is 
introduced when solving the Schrodinger Equation. 

For the 7=1 5=2 case of B*B* system, we found, in 
all the three SU(3) quark models, a exchange provides 
attraction but a' and tt exchanges provide repulsion, and 
the total attractive interaction is too weak to bind the 
B*B*. 



C. DD* 

The only difference between the systems of DD* and 
BB* is the flavor of heavy quark . Considering the resem- 
blance between DD* and BB*, the one meson exchange 
interaction potentials should have similar properties. In 
our calculation we find the one meson exchange inter- 
action potentials of DD* system are a bit larger than 
those of BB* system with the same quantum numbers of 
isospin 7 and C-parity C, and the total interaction also 
is a bit larger. That is to say, the heavy quark mass has 
only a little effect on the interaction potentials (see Eq. 
19). After solving the Schrodinger Equation, we list our 
resuhs in Table [TTT1 

From Table IIIIl we find the result of DD* system is 
quite different from that of BB* system. This is because 
the much lighter reduced mass makes the DD* system 
more difficult to form a bound state. We see that only one 
5- wave 0"''(1"'"'") DD* bound state with binding energy 
of 12.72-47.39 MeV could exist. Unlike the BB* system. 
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TABLE III: Binding energies (MeV) of the bound states of 
DD* system with different / and C quantum numbers. Sym- 
bols have the same meaning as in Table H] 



/ 


c 


I 


II 


III 





+ 


12.72 


47.39 
24.26 


33.32 
14.24 


1 


+ 









the total attractive interactions are too weak to form 
0-(l+-) and DD* states. 

In Fig. 13 the total interaction of /=0 C=+ DD* sys- 
tem is shown. Our 0"'"(1'^'*") DD* bound state has a mass 
3832-3867 MeV, and it fits well with the mass of new res- 
onance of X(3872) observed by (26l - [29| . Particularly, the 
quantum numbers of 0"'"(1"'""'") seem to be more prefer- 
able in the recent experiment [i^ of X(3872). Here we 
speculate that the study of the 0+(l+"'") DD* molecule 
bound state may shed light on the X(3872) structure. 




-200 - 
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-300 I ' 1 ' 1 ' 1 ' 1 ' 
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FIG. 3: The total interaction potential of I=Q C=+ DD 
system. Symbols have the same meaning as in Fig. [T] 



D. D*D* 

In D*D* system the one meson exchange interactions 
have similar properties to those of B*B* system. Analo- 
gous to the process of B* B* system, we list the result of 
D*D* system in Table |lVl 

Comparing Table |IV] with Table HH we find the D*D* 
system is more difficult to form bound states than the 
B*B* system. This is because the reduced mass of the 
D*D* system is 1003 MeV which is lighter than that of 
B*B* system of 2663 MeV. For D*D* system, we see 
there is only one 0+(2++) D*D* bound state with bind- 
ing energy of 14.33-50.22 MeV. Unlike the B*B* system, 
the total attractive interactions are too weak to form 
0-(l+-), 1"(0++) and 1+(1+-) D*D* states. 



TABLE IV: Binding energies (MeV) of the bound states of 
D* D* system with different isospin / and spin S. Symbols 
have the same meaning as in Table IH 



/ 


s 


I 


II 


III 









16.15 


7.97 





1 




26.38 


15.89 




2 


14.33 


50.22 


35.69 


1 



1 
2 









The total interaction potential of /=0 S^2 D*D* sys- 
tem is depicted in Fig. H This 0+(2++) D*D* state 
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FIG. 4: The total interaction potential of /=0 S=2 D* D* 
system. Symbols have the same meaning as in Fig. [T] 

has a mass of 3964-3999 MeV. We notice that the mass 
and quantum numbers of this state are consistent with 
those of y(3940) observed recently by Refs. [13^1 and 
discussed by Refs. [U SI]. Therefore, the study of our 
0+(2++) D*D* molecule state might be helpful to un- 
derstand the structure of the y(3940). 



IV. SUMMARY 

In this work, we have performed a systematic in- 
vestigation of the bound state problem of DD*,D*D*, 
BB* and B*B* systems. In DD* system, we find only 
one O"*" (!"'""'") bound state and this molecule state might 
partly account for the structure of Ar(3872). In D*D* 
system, we find only one 0+(2++) bound state could ex- 
ist and this D*D* molecule state might be explained as 
V(3940) observed in experiments. 

General speaking, we found that BB* and B*B* can 
form bound states more easily than DD* and D*D*, be- 
cause of their much heavier reduced masses. For the two 
charged new resonances of 2';,(10610) and 2';,(10650), our 
calculation shows that they might be explained as 1+(1+) 
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weakly bound states of BB* and B*B* , respectively. For 
BB* system, other two bound states are found in our ap- 
proach as 0+(l++) with a mass 10506-10557 Me_V and 
0-(l+-) with a mass 10538-10600 MeV. For B*B* sys- 
tem, other three bound states are also found. They are 
0+(2++) with a mass 10552-10603 MeV, 0^(1+- ) with 
a mass 10583-10645 MeV, and l+(0+-) with a mass 
10646-10649 MeV. However, until now no experimental 
evidence has been reported for these five bound states in 
our approach. 

What's more, there are some states unbound in our 
chiral SU(3) quark model (model I) but bound in our 
extended chiral SU(3) quark model (model II and model 
III). So far, we can't draw definite conclusions for them. 



It is expected future experimental measurements will pro- 
vide more information about these new hadron structures 
and provide a discrimination for our model calculations. 
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